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Abstract 

A new method of modeling of the current signal induced by charged 
particle in silicon detectors is presented. The approach is based on 
the Ramo-Shockley theorem for which the charge carrier velocities 
are determined by taking into account not only the external electric 
field generated by the electrodes, but also the Coulomb interaction 
between the electron and hole clouds as well as their diffusion. 
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1 Introduction 

It is obvious that identification of particles and fragments produced in nuclear 
reactions is crucial for any experimental or technical work in nuclear physics. 
Among different ways of identifying charged particles the classical AE — E 
telescope method remains still the flagship. Recently, an alternative method 
based on the Pulse Shape Discrimination (PSD) technique applied for silicon 



detectors is being developed and is increasingly drawing attention. Recent 
results demonstrate that the method can offer charge and isotopic identifi- 
cation comparable to that obtained with the classical AE — E method. The 
main advantage of the PSD method comes from the fact that it requires 
only one electronic channel for detection and identification. It is thus an 
important point for designing and constructing multi-detector systems. 

A significant difference between the AE — E and the PSD techniques 
results from the fact that the former is governed basically by the energy loss 
process (Bragg curve), while the latter is primarily related to the Plasma 
Delay Effect (PDE) [1-6]. In silicon detectors, this effect manifests itself with 
shortening of the pulse rise time with decreasing Z for low and intermediate 
mass fragments, for which the generated charge is practically completely 
collected by the detector electrodes. The experimental data demonstrates 
that the PDE concerns particles with small Z, for which the Pulse Hight 
Defect (PHD) is still of little importance. 

For better understanding of the identification idea associated with PSD 
technique, and for its future development, it is crucial to have at ones com- 
mand a perfect simulation of the time dependence of the experimental signal 
produced by an ion with a given charge, Z, atomic mass, A, and energy, E. 
The main goal of such a simulation is to describe the extraction and collection 
of the generated charge carriers moving in the external electric field distorted 
due to the presence of a highly ionized track and due to the diffusion process 
of the carriers. 

As presented in [7], an approach in which the distortion of the electric 
field caused by the generated carriers is neglected, is able to correctly describe 
the current signals for light charged particles (LCP), e.g. protons. However, 
this simplified approach completely fails in case of heavy ions (HI) for which 
the collection time of the generated carriers gets longer (t hi > t lcp ). His- 
torically, this difference, associated to a slower carrier collection for HI, was 
quantified as a plasma delay (PD) effect. Since this effect influences the 
current signal rise-time, it appears to be crucial for the PSD technique. 

An attempt to describe phenomenologically the delayed carrier collection 
time in silicon detectors has been recently proposed in [8]. The proposed 
description took into account the polarization of the electron-hole pairs gen- 
erated by the HI and connected it to the relative dielectric permittivity. 
Another important assumption was that the dissociation of pairs in time oc- 
curred with a constant probability and the modified electric field, inside and 
outside of the ion range, was given by the Maxwell equation for the electric 
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field in the inhomogeneous medium. With these assumptions, the model was 
indeed able to describe the experimental pulse shapes quite accurately. 

In the present paper we propose another, more microscopic approach. 
The main model assumptions are the following: 

i. Propagation of the electric charges (electrons, holes) gener- 
ated in the detector is represented by evolution of the Gaussian 
clouds for which the centroids and variances are treated as inde- 
pendent variables. 

ii. Position of the centroid of each Gaussian is governed by 
the drift process, while its variance undergoes both, diffusion and 
the drift process. 

The first results of the model calculations indicate some binding effects be- 
tween the holes and electrons, in a region similar to that predicted by the 
phenomenological model of Ref. [SJ. In this region also the electric field 
shows similar behavior to that presented in [8j. 

The detailed description of the new model is presented in the following 
section. The preliminary results of the calculations and comparison to the 
experimental data are presented in Section 3; Conclusions and possible ex- 
tensions of the model applicability are given in Section 4. 



2 Description of the model 

A particle entering the silicon detector is assumed to degrade its energy 
according to the Bragg curve which relates the generated ionization B(x) to 
the particle position x. We assume that the X direction is perpendicular to 
the detector surface. In order to describe the initial, local, density of the 
electrons p e {v,t = 0) and holes Ph(r,t = 0) we assume that the ionization is 
proportional to the local stopping power B{x) = — ^-(x), where w = 3.62 
eV is the energy for an electron-hole pair production and jr-(x) is the local 
stopping power [9\. Just after stopping of the impinging ion, the carrier 
density can be described as: 



p e {r,t = 0) = -p h {r,t = 0) = - J B{x')6{x - x')5{y)5{z) dx' (1) 

where r = [x, y, z\. This assumption states that for t — the ionization is 
localized along the X axis only and disappears elsewhere. 



3 



In order to describe the time evolution of the generated ionization we 
assume its distribution in the following form: 

p e (r,t) = J B e (x',t)G e (x-x',y,z,t)dx' (2) 



p h (r,t) = j B h (x',t)G h (x -x',y,z,t)dx' (3) 
which is analog to ([I]) and we set: 

- B e (x, t = 0) = B h (x, t = 0) = B(x) (4) 
Functions G e and Gh are assumed to be Gaussians: 

n ( ' a 1 ( (x - xf + y 2 + z 2 . 
G e (x-x ,y,z,t) = —== exp — — (5) 



and 



G h (x - x',y,z,t) = — 1- exp I — — I (6) 



If a e — > and ah — > for t — > then the functions G e and Gh can be 
regarded as representations of the 5 function, thus: 

G e (x — x', y, z, t — 0) = Gh(x — x', y, z, t — 0) = 6(x — x')8(y)5(z) dx (7) 

Now, the goal is to describe the time evolution of the functions B e , Bh and 
G e , G e (for determination of G e , G e it is sufficient to derive the time evolution 
of their variances a\ and a 2 ,). In order to do it we define the one dimensional 
densities associated with the variable 



oo oo 



r) e (x,t) = J dy J dzp e (r,t) (8) 



- oo — oo 



OO OO 



r) h (x,t) = J dy j dzp h (r,t) (9) 



- oo — oo 
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and we divide the thickness of the detector dsi into N intervals Ax = 
Let us assume that in the interval Axi (i — 1, N) the associated values of a e i 
and Ohi do not change substantially within the radius of a few sigma around 
Axi, and that a linear approximation can be used for the functions B e (x',t) 
and B h (x',t) within Axf 

B e (x',t)=p ei (t)x' + q ei (t) (10) 

B h (x',t)=p hi (t)x' + q hi (t) (11) 

With the above assumption for x within an interval Axi one can approximate 
the densities: 



oo oo 



r} e (x,t) ~ J dx j dy J dz (p ei (t)x' + q ei (t)) G e (x-x' ,y, z,t) = p ei (t)x+q ei (t) 

oc 

(12) 



-oo — oc 



OO OO OO 

r] h (x,t)^ J dx' J dy J dz (p ei (t)x' + q ei (t)) G h (x-x' ,y, z,t) = p hi (t)x+q hi (t) 

— oo — oo — oo 

(13) 

which means that, in practice, one can use the same coefficients for linear 
expansion of both, the densities r) e , r\ h and of the functions B e , B h . 

We introduce also the x e0 (t) and x e jy(t) coordinates, which denote the 
beginning and end of the B e (x,t) distributions for electrons. Similar coor- 
dinates Xho(t) and XhNif) are introduced for holes (see right-upper panel on 
Fig. 2). 



2.1 Electric field determination 

In order to determine the drift velocity associated with the centers of Gaus- 
sians G e , Gh we have to calculate the respective effective electric field acting 
on the carriers described by above distributions. Such a field is determined 
by a static voltage applied to the detector electrodes and by the Coulomb 
interaction between the Gaussian charge clouds. The detector static field at 
position x, considered from the rear side of the detector (order n — p from 
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the point of view of the particle entering the detector, see e.g. 0), is given 
as: 



2V d x , V-V d .... 

E st at{x) = -jj- + — (14) 

a Si d Si 

where the bias voltage V is assumed to be higher than the depletion voltage 
Vd which for the bulk concentration of donors Np and permittivity 
reads as: 

V, - (15) 

In order to find the modification of the electric field caused by the gen- 
erated plasma, let us consider two Gaussians describing the distribution of 
the charges Zi, Z 2 , centered at a relative distance r 12 =| r i ~~ r 2| an d char- 
acterized by variances af, a\ respectively. The mutual interaction potential 
of the clouds can then be expressed in the form which can be often found in 
quantum molecular dynamics calculations (see e.g. [TO] ) 

(ri-rpi) 2 \ / — frg— £0g) 



, e 2 z x z 2 rr ex P { 2a t ' ex P , , 3 

< r oi, ^02, o u cr 2 ) = — -3 // : : d r x d r 2 

(277(7! a 2 ) JJ ki-r 2 | 



,2 



erf (7i 

Z,Z 2 V ^°^ (16) 

where a = \fo\ + o\. 

Let us now assume that the intervals Axi are small enough to enable the 
linear approximation for the charge densities p e , ph and the functions B e ,Bh 
with the use of the coefficients p and q. For simplicity, we introduce variables 
p , q in lieu of the p e i(t), Puif) and q e i(t), qmif)- If we denote the endpoints 
of the Axi interval by c and d then for a Gaussian centered at a point a 
and representing the charge Z a , its interaction with the charge located in the 



interval Axi = (c, d) given by (12) and (13), can be formulated as: 

e 2 f d erf (t!tO 

V c {a,c,d,p,q,a a ,ai) = —Z a dx (px + q) — ^— (17) 

s J c \x — a\ 

where the a a and <7j above denote the standard deviations of the Gaussian 
describing the charge Z a and of the Gaussian from the interval Axi, respec- 
tively, and cr s = \J o" 2 -4- of. 
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The above form allows us to describe the respective effective electric field 
E x acting on the Gaussian located at a point a as 

\x—a\ 



1 OV e 2 f d d er ^ ( ) 

E x (a,c,d,p,q,a a ,<Ji) = - — = / dx (px + q) - . ^— (18) 

Z a oa e J c oa \x — a\ 

Since J~ = — the above formula can be expressed as: 

8 «■* 



e 2 r 

S x (a, c, d, p, g, a a ,ai) = — dx (px + q) 

£ J c 



dx \x — a 



2 / erf ( ^1 1 erf I I \ 2 /-d a erf I 1^ 



q —r, f ; s r— H p dxx- . s — (19) 

e I a — a c — a / £ ./„ arc kr — a 



c 



After integrating by parts one obtains: 



2 / erf ( J^l ] er f f J£zM \ \ 2 / erf ( I erf I ^1 



£ \ a — a c — a / £ \ a — a c — a 



,2 r d erf 



/ |g-a| 



--p ( d dx ~\ KV ~ 2 y (20) 

£ y c |x — a\ 

The last integral can be easily evaluated by expanding the error function. 

The total effective electric field modified due to the presence of plasma can 
be obtained by summing up contributions of charges located in all intervals 
(cj, di) and of their mirror charges induced in the detector electrodes (see 
Fig. 1). 

In order to determine the time evolution of the charge distribution one has 
to know, in addition, the generalized force associated with the a a variable. 
Similarly as for the effective electric field, the net value of the force acting in 
the a a direction can be obtained by summing over all ingredients associated 
with the charge distribution. This leads to the following expression for the 
interaction with the charge located in the Axi interval: 

9Vc _ e 2 v , ^ d erf (S) _ 



dV c e 2 f a 

7 a a ,a,cd = = ~~ Za J dx (P X + q) 



da a \x 
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= ~ a 7%^j I exp I I ~ exp 



V 2 

+ Z „^±i)(e rf (^))-e rf (^))) ( 21 , 

Now we are in position to calculate the time evolution of the charge 
generated in the detector. This process is determined by the drift and by the 
diffusion of the interacting clouds of electrons and holes. We will consider 
the evolution of the centroids and of the variances of Gaussians representing 
a fraction of the charge distribution located in the middle of the intervals 
Axi and at the start- and end-points of the distributions of electrons and 
holes (points x e0 (t), x eN (t) and x h0 (t), x hN (t)) . 

2.2 Evolution of the function B 

In the present subsection we will describe the numerical method used to 
determine the time evolution of the ionization clouds. In the following we 
assume that the evolution of the functions B e , Bh is determined by the time 
evolutions of the coefficients p ei , q ei and p^, qu- In order to find how the 
expansion coefficients p ei , q ei and pu , qu propagate in time, we have to 
investigate the time evolution of the functions i] e (x,t) and i] h (x,t). 

Below we consider the formulas for electrons only, keeping in mind that 
the formulas for holes are analogical. 

As we will show later 

d d 

— r} e (x, t) = -— (r} e (x, t)v xe (x, t)) (22) 
thus, the differential dr] e (x,t) can be written as 

( d d \ 

dr) e (x,t) = - lv xei (x,t)—r) e (x,t) + r) e (x, t)— v xe (x, t) 1 dt (23) 

If the Xi denotes the center of the interval Ax^ and q e i and p e i are the 
coefficients of linear expansion 

r} e (x,t) = Pei(t)(x - Xi) + q ei (24) 
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and, if one denotes the average velocity and the average linear density asso- 
ciated with the interval Axi by v xe i(t) and f] xe i{t), respectively, then: 

dijei(t) = - (v xei (t)p ei (t) + r} ei (t)(p xei (t)) dt (25) 

Here, ip xe i(t) is the differential coefficient of v xei (x,t) at the point Xj. As one 
can see, in order to calculate the above increment, we have to trace the time 
dependent tables of rj e i(t), v xe i(t). The tables of p e i(t), q e i(t) and <p xe i(t) are 
obtained by fitting the smooth curves to the distributions rj e i(t), v xe i(t) in 
every time step. For t — the rj e i(t = 0) is given by the Bragg curve. In 



order to make use of the formula (25) we need to construct the respective 



tables for velocities. Knowing the effective electric field for electrons in the 
x direction, E xei , one can assume that the respective average velocity of the 
center of Gaussian located at a point Xi is proportional to the strength of the 
field: 

V X ei f^xe-^xei (26) 

where [i xe and fi x h are the electron and hole mobilities, respectively. 

Knowing the drift velocity, one can calculate the evolution of the charge 
deposited in every interval Axi including the edge intervals with variable 
ends x e0 (t), x eN (t) and x h0 (t), x hN (t). 

2.3 Charge propagation in the perpendicular direction 

The diffusion and transport processes in the electric field influence also the 
widths of the charge distributions located in every Axi interval. Extending 
the above reasoning we can assume that the velocity, v aei , describing the 
rate of the standard deviation expansion in perpendicular direction has three 
components: 

Vaei^V^ + V^ + V^ (27) 



The first term, v„ ei , results from the field described by (21). In analogy 
to the charge drift in x direction, one can assume that this component is 
proportional to the field acting on the charge Z a associated with the Gaussian 
with a standard deviation a a 



E ae i = (28) 
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where F aa is the net force given by interaction (21). Thus, the respective 
velocity can be expressed as 



v E aei = ^ ae E aei (29) 

where, the \i ae parameter is the only free parameter of the model. It seems, 
however, that it can be determined theoretically in the future. A similar 
parameter for description of the hole propagation can be calculated assuming 
the following proportion: 

Mae = M*e ^ 

The velocity v^ ei follows from the solution of the second Fick's law for diffu- 
sions of Gaussian density distributions: 

d da ei D e 

V ™ ~ ~df ^=0,v- ei =0- ^ ( 31 ) 

where D e is the diffusion coefficient for electrons. 

Another process which affects, on average, the widths of the respective 
Gaussian distributions used to describe the charge located in the interval Axi, 
is related to the transport of the carriers. In order to describe this process 
we consider an increase of the variance, a 2 , of the Gaussian distribution of 
the charge of Z particles contained in an interval Ax, with linear density 
rj(x) = Let the ^r 2 denote the sum of squares of deviations of particle 
positions from the average. In this consideration we neglect the influence of 
the diffusion and of mutual interactions of clouds on the propagation of the 
variance. 

v 2 

The change of a 2 = , resulting from the flow of particles into and 
out of the cell Ax (with the net value of dZ) is, in general, equal to 

, 2 d(£rf) a 3 (x)dZ , , 

da 2 = -^f^- - -^j— (32) 

As one can see, in order to find the increment da 2 we have to find the 
increments dZ and d (Y^ rf) . 

Let us begin with the description of d rf) . If the accretion of particles 
in an interval As across the points x\ = x — Ax/2 and X2 = x + Ax/2 is 
denoted by d^and dZ 2 , respectively, and the variances at these points are 
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denoted by o\ = o 2 {x — Ax/2) and o\ = a 2 (x + Ax/2), respectively, then 
the increment of the sum r 2 can be determined as: 



«(E 



a\dZ\ — o\dZ\ 



(33) 



Denoting the velocities of particles at points X\ and x 2 by v(x — Ax/2) and 
v(x + Ax/2), respectively, the accretions dZ\ and dZ 2 can be determined as 



dZ 1 = rj(x - Ax/ 2) v(x- Ax/2)dt 



(34) 



dZ 2 = 7](x + Ax/2) v(x + Ax/2)dt 
Now one can calculate the increment d rj) as 



(35) 



4E 



-dt ■ Ax- 



a 2 (x + Ax/2)rj(x + Ax/2)v(x + Ax/2) — o~ 2 (x — Ax/2)?](x — Ax/2)v(x — Ax/2) 

Ax 

(36) 

The expression in square brackets tends to the partial derivate d ^ a ( X W X ^ V W) 
for Ax 0. 

Similarly, the increment dZ = dZ\ — dZ 2 can be written as 



r)(x + Ax/2)v(x + Ax/2) - n(x - Ax/2)v(x - Ax/2)] , 

dZ = —Ax dt 

Ax 

(37) 

and again, in the limit of Ax — > the expression in square brackets ap- 
proaches to dO?(zMz)) ^ -j^Qg already been used in (22). 



Taking the above into account, setting Z = r/(x)Ax and taking the Ax — > 



limit, eq. (32) can be transformed into: 



da 2 (x) 



—d (a 2 (x)rj(x)v(x)) a 2 (x)d(r](x)v(x)) 



dx 



+ 



dx 



dt —v(x)d (a 2 (x)) 



T)(x) 



dx 



dt 



what gives 



(38) 
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da 2 da 2 da da , . 

=► »r = (39) 



<9t <9x <9t <9x 

Finally for = and v® ei = one can write 



V ™ = ~7u~ \v E =0,v D =0= -Vxei^- (40) 

at aes ox 



In order to use the above formula we trace the changes of the vector ofcr e 

E . V D . 



values as a function of the position index, i. Knowledge of velocities v^ ei , v D 



and uj ei allows to calculate the propagation of the width of the distribution 
of electrons. The formulas for holes are analogical. 



3 First prediction of the model and compar- 
ison with the experimental data 

For the first comparison of the model prediction with the experimental data 
we choose the data for l2 C ion which have already been used in [8]. This 
gives also the opportunity to compare the present model predictions with 
those obtained in a more phenomenological approach. For the measurement 
the neutron transmutation doped (n-TD) silicon detector [TT] was used. This 
n-type bulk and extremely thin p-type zone has a thickness dsi = 310 fim. 
The energy measurement was performed using the charge output, while the 
current pulses were measured using the current output of the same pream- 
plifier, described in [7]. This paper presents also in detail the experimental 
setup and conditions used for the 12 C ions (and LCP). 

Before describing the induced current pulse, we focus first on the propaga- 
tion of the electric field and the propagation of the electron and hole densities 
in parallel and perpendicular directions. The evolution of these observables 
is important for understanding the mechanism of the plasma delay process. 
In the following, we consider an 12 C ion impinging on the n-type rear side 
of the silicon detector. This, so called "rear-mount", gives quite different 
shapes as compared to the "standard mount", and these pulse shapes are 
much better suited for the PSD technique [7]. 

For actual calculations it is necessary to set some physical coefficients 
describing the electric field propagation, as well as coefficients describing the 
drift and diffusion process in silicon. In table 1 we collect values of these 
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parameters: 



feature 


symbol 


value 


remarks 


energy per e-h pair creation 


w 


3.6 


eV 
pair 




material constant 


silicon dielectric, permittivity 


e r 


1] 


..7 


material constant 


electrons mobility 


f-^xe 


135 


am 2 
Vns 




material constant 


holes mobility 




47.5 


/im 2 
Vns 




material constant 


electrons variance mobility 




r> am 2 

Vns 




free parameter 


holes variance mobility 


Hah 


f^ah Mere ,, 


model assumption 


diffusion coefficient for electrons 


D e 


3.49 


am 2 
ns 




material constant 


diffusion coefficient for holes 


D h 


1.228 


um 2 
ns 




material constant 



As already mentioned, the electric field propagation results from the static 
detector bias and from the generated charge density propagation. At the 
starting point, when the generated electrons and holes are almost exactly at 
the same positions, the electric field is still equal to the external one (14). 
This field, for t=0, is denoted in Fig. 2 by a doted line. After this initial 
moment the static field causes the shift of the electron and hole distributions 
and therefore in the next moment some of the carriers are moved outside 
of the overlap region. Next, the variance of the Gaussian partial density 
distribution associated with these carriers begins to grow, due to the non- 
compensated electric field in perpendicular direction. This effect, which is 
displayed in Fig. 3, causes breaking of the initial bonds between the electrons 
and holes. As a result, the considered electrons and holes start leaking slowly 
from the overlap region and begin moving in opposite directions. At the same 
time, the increasing shift between electrons and holes leads to significant 
reduction of the electric field in the interaction region. This scenario is well 
associated with the postulates of the phenomenological model of ref. [8]. 

The evolution of the charge density, the effective electric field and the 
variances of the Gaussian partial densities are presented in Figs 2-5 in which 
the blue lines represent the dependences for the electrons while the red ones 
represent the holes. Fig. 2 shows the evolution of the linear electron and hole 
densities. One can see that, during the first 20-30 ns, the electrons and holes 
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remain bound in the region close to the detector surface. After that time, 
one can observe the electric field restitution practically in the whole detector 
area. The time behavior of the electric field is presented in Fig. 3. In order 
to save the calculation time, the field is calculated only at points where the 
density of the particles is not equal to zero. Fig. 4 presents the evolution 
of the width of the charge distributions. We can notice that in the overlap 
region the effective electric field is reduced to very low value. As one can see, 
the evolutions of the electric field and of the charge densities (in parallel and 
perpendicular directions) are strongly correlated. 

Up to now, the mutual Coulomb interactions between the charge clouds 
have been taken into account. In order to see the importance of these mu- 
tual interactions, they have been neglected in the charge density evolutions 
presented in Fig 5. As one can see, in this case, the collection time becomes 
about tree times shorter, due to the lack of binding between the electrons 
and holes. Knowledge of the charge propagation, by using of the Ramo- 
Shockley theorem [T21 [321 EH] , allowed as to determine the current pulse time 
dependence. 

In order to obtain a rough estimate of the pulse shape, we describe the 
partial current associated with the Gaussian cloud in approximate way. For 
simplification, using Ramo-Shockley theorem, we replace the Gaussian charge 
distribution by a point-like one. Such an approach neglects effects associated 
with the charge diffuseness, particularly for clouds moving closely to the 
detector electrodes. 

The result of such a calculation is presented in fig 6. The total current 
pulse is denoted by black solid line while the electron and hole contributions 
are represented by the blue and red ones, respectively. 

For comparison of the calculated pulse with the experimental one the 
primary pulses from Fig. 6. have been corrected (see [3]) for the preampli- 
fier's response. The results are shown in Fig. 7. We have to stress that in 
the present calculations we did not search for the best value of the fi ae pa- 
rameter. We also did not consider some quite complicated factors, specified 
below, which could affect the obtained results and which will be a subject of 
the forthcoming paper: 

i) precision of the Energy-Range tables (average accuracy of about 10%, 

ii) diffuseness of the Gaussian clouds and its presence in the application 
of the Ramo-Shockley theorem, 

iii) dead layers of the detector and their effect on the measured energy (as 
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one can see on Fig. 8, for 80 MeV 12C the collection time is very sensitive 
on the ion energy), 

iv) accuracy of the active detector thickness and of the electric field de- 
termination, 

v) accuracy of the preamplifier response description. 

In order to demonstrate that the present model is able to describe cor- 
rectly the plasma delay effect, in Fig. 8 we present the correlation between 
the collection time and the energy loss for 10 B , 12 C , and 14 N ions. We use 
the collection time rather than the experimentally preferred rise time, noting 
that these two observables are strongly correlated. Fig. 8 shows that the 
model can reproduce the experimental trends, especially the characteristic 
"back-bending" of identification curves at low energies. 

4 Conclusions 

We have proposed a description of the evolution of charge density and of the 
effective electric field by taking into account the mutual Coulomb interactions 
between the charge carrier distributions. According to the present approach 
the plasma delay effect is associated with the propagation of the carriers in 
both directions, perpendicular and parallel to the primary ionization path. 

The duration of the obtained pulse for the 80 MeV 12 C ion, corresponds 
quite well to the one obtained in the physical measurement. Also the shape 
of the Energy-Collection Time correlation and its element dependence are 
quite well reproduced by the model. Nevertheless, the model still needs to 
be confronted with a broader collection of the experimental data obtained 
for detectors of various thicknesses and biased by various voltages. 

Once tested on a broader collection of the experimental pulse shapes, the 
model will enable the theoretical search for the best identification method 
based on the pulse shape analysis. It will also enable the study of the depen- 
dence of the identification quality on the detector thickness and maybe on 
some special construction of the detector with non-linear electric field (ob- 
tained by the heterogeneity of doping) for regions with poor resolution (small 
ion energy, see Fig. 8). The presented approach should be also suitable for 
testing the temperature dependence (via the respective dependence of the 
diffusion coefficients). 

In order to draw some more quantitative conclusions from the comparison 
of the model with the experimental data one has to estimate the uncertainties 
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related to various possible ingredients, mentioned in the previous paragraph. 
Also we have to make an attempt to constrain the single free parameter \x ae 
from the respective classical consideration. 

The actual model calculation is quite time consuming. For standard 
processor the calculation of one pulse associated with the 80 MeV 12 C ion, 
takes about 5 hours of CPU, thus some code optimization is still needed. 
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Figures 

Figure 1. Allowing the charge induced by plasma in the detector electrodes. 
In the present approximation we consider only the nearest mirror reflections. 
So the influence of the charge induced in an electrode positioned at on the 
detector electric field, acts as a part of the Gaussian localized in position — x 
with the reverse charge and with the same variance. Similarly, the charge 
induced in the second electrode is represented by a respective Gaussian lo- 
calized in a symmetric point at 2dsi — x . 
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Figure 2. Propagation of the linear density of the electrons rj e (blue lines) 
and the holes r]h (read lines) due to the ionization induced by 80 MeV 12 C 
ion entering the Si detector from the rear side. 
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Figure 3. The effective electric field strength inside the silicon detector 
at different moments in time, due to the ionization induced by 80 MeV 12 C 
ion penetrating the detector from the rear side. The doted line gives the 
undisturbed electric field for t = 0. The field is calculated only at points 
with coordinate x where the the density of electrons (blue line) and holes 
(red line) is not equal to zero. 
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Figure 4. Time evolution of the width of the charge distribution which 
determines the charge propagation in the perpendicular direction. Standard 
deviations for electrons are represented by the blue line and variances for the 
holes by the red one. 
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Figure 5. Same as Fig. 2, but neglecting the mutual Coulomb interactions 
between the carrier clouds. Significant difference in the charge collection time 
can be observed as compared to the complete case (see Fig. 2). 
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Figure 6. Model prediction of the current signal induced by an 80 MeV 
12 C ion penetrating the silicon detector from the rear side. The predictions 
are not corrected for the preamplifier response. The solid, black line repre- 
sents the total signal, the blue line presents the electron contribution while 
the red one the hole contribution. Mean experimental current signal is pre- 
sented by the dashed line. 
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Figure 7. Same as Fig. 6 but the respective lines have been corrected for 
the preamplifier's response. 
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Figure 8. Model prediction for correlations: Energy vs Collection Time. 
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